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Abstract 

On 

Recently in boson star models in framework of Brans-Dicke theory, 

three possible definitions of mass have been identified, all identical in gen- 
eral relativity, but different in scalar-tensor theories of gravity.lt has been 
conjectured that it's the tensor mass which peaks, as a function of the 
central density, at the same location where the particle number takes its 
maximum. This is a very important property which is crucial for stability 
KS ' analysis via catastrophe theory. This conjecture has received some nu- 

C^ ■ merical support. Here we give an analytical proof of the conjecture in 

framework of the generalized scalar-tensor theory of gravity, confirming in 
this way the numerical calculations. 

PACS numbers: 0440D, 0450 



Boson stars were first discussed by Kaup |]T] and then by Ruffini and Bonazzola 
0. Boson stars in scalar-tensor theories of gravity have been investigated exten- 
sively by many researchers. The first model of a boson star in pure Brans-Dicke 
theory has been studied by Gunderson and Jensen ||. Their work was gener- 
alized by Torres |3| who was studied boson stars in scalar-tensor theories with 
non-constant u>bd{&) -More recently boson stars have investigated in the papers 
by Torres et al. ||,f| and in the paper by Comer and Shinkai 0. In || boson 
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stars have been studied in connection with so-called gravitational memory ||, 
while their stability through cosmic history has examined in |7| and M. Finally, 
the dynamical evolution of boson stars has been investigated in the paper by 
Balakrishna and Shinkai [f|. For more details we refer the reader to the most 



recent review on boson stars [ID]. 



Here we consider complex scalar field boson stars in the most general scalar 
tensor theory of gravity with an action in Jordan frame 



S = --tV / \f-$ i F ^)R - H($)r u d^d v $ + &((/>)) d 4 x + (1) 

107TG* J V V ' 

+ f \/-£ (-zg^d^+d^ - W(tf + tf )) d*x 



where R is Ricci scalar curvature with respect to the space-time metric g^ U) G* is 
the bare Newtonian constant and $ is the gravitational Brans-Dicke scalar with 
potential term U((f>). \P is massive self-interacting complex scalar field with 

Hereafter we will consider only static and spherically symmetric boson stars. 

The explicit form of the action shows that we have [/(l)-invariant action 
under global gauge transformation \1/ — > e m \l/, a being a constant. This global 
L r (l)-symmetry gives rise to the following conserved current 

^ = \4- Z 9~9 IMV (^d^ + - V + d^) (2) 

The conserved current leads to a conserved charge - total particle number making 
up the star 

N= f J°d 3 x (3) 

Binding energy is then defined by 

E B = Mstar - mN (4) 

where m is the particle mass. 

Now a problem arises: How to define the mass which appears in the expression 
for the binding energy? 

In contrary to general relativity the definition of mass in scalar-tensor theories 
of gravity is quite subtle. This problem has been recently examined numerically 



in the work by Whinnett ITTJ.He has considered three possible mass definitions in 
Jordan frame, namely Schwarzschild mass Ms (i.e. ADM mass in Jordan frame), 
Keplerian mass M^epier and the tensor one My. As it has been shown numeri- 



cally in |ITT[] (see also [13]) these three masses differ significantly from each other 



in the case ujbd = — l.Keplerian mass leads to positive binding energy which 
means that every boson star solution is in general potentially unstable. Contrary, 
Schwarzschild mass leads to negative binding energy suggesting that every solu- 
tion is potentially stable even for large central densities p = | ^(O) | 2 . 

It should be noted that for large constant ujbd (say ujbd > 500) the difference 
between the three masses is negligible. However, for arbitrary 1 uobd{4>) it's possible 
that the three masses may differ from each other significantly. This may occur in 
the early universe when the cosmological value $oo is sufficiently smaller than 1 
||. Moreover, our numerical calculations show that for some physically relevant 
functions ubd^) we may have My — Ms ~ (0.15 — 0.20)My. On the other hand 
the scalar tensor theories of gravity with ujbd — ~ 1 better describe the early 
universe(see |13| and references therein). So, the case when ujbd is not large, is 
also physically relevant. Therefore, when we study the boson stars in the early 
universe the mass choice is crucial. 

It's the tensor mass which leads to physically acceptable picture. In jnj it's 
shown numerically that the tensor mass peaks at the same point as particle 



numbers - a very important property in general relativity [14]. This property is 
also crucial for the application of catastrophe theory to analyze the stability of 
the boson stars |L5| , |Tj| . 



In [IT | and [0J the problem for analytical proof of the fact that it's the tensor 
mass which peaks at the same location as the particle number has been stated. 
The main purpose of this paper is to fill this gap. 

In our opinion it's more convenient to work in Einstein frame given by 

g»» = F($)gp, (5) 

In Einstein frame the action ([|) takes the form 

s = -— Jg- JV=g(R- 2<r^<R0 + u(4>)) d 4 x + (6) 

where R is the Ricci scalar curvature with respect to the metric g^, U{(f>) = 

A 4 (0)[/($(0)) and A 2 {(t>) = F" 1 ($((/))) as <j> is given by 



di 
The action (0) leads to the following fields equations 



3f dln(F(Q) V lff(Q 
N4 V ^ J + 2F(0 (7) 



g; = k*t; + 2d^<f> - d^d^s; + ~u(<f>)6; 



x We mean arbitrary ojbd(4>) f° r which the theory passes through all known gravitational 
experiments. 



D0+I/7'(0) = -^a(0)T 



□tf + 2a(<p)d a (f)d (7 y = -2A\cj)) 



, ^dW(ty+ty) 



m+ + 2a{<f))d a (j)d (T ty + = -2A 2 (<j)) 



., L .dW(V + V) 



where □ is d'Alambert operator in terms of the metric g^ u , k* = 87rG*, «(</>) = 
-jT ln(A((j))) and T is the trace of the Einstein frame energy- momentum tensor of 
the complex scalar field given by 

t; = ^a 2 ^) (d^ + d^ + d^d^ + ) - (9) 

--A 2 (<P) (d^+d^ - 2A 2 (<j))W^ + ^)) 5; 

The conserved t/(l)-current in Einstein frame is 

,p = -A 2 (<p)^g» u (ttd„*+ - V + d„V) (10) 

As we have already mentioned we consider static and spherically symmetric 
boson stars i.e. space-time with a line element in Einstein frame 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 (0)tfy> 2 ) (11) 

and complex scalar field in the form 

m = a(r)e iujt (12) 

where u is real positive number and a(r) is a real function. It should be noted 
that the asymptotic behaviour of the field (p set the following tight constraint 

< u < m (13) 

In this case the field equation system reduces to the following system of ordinary 
differential equations 

A' = ^—^- + /€,e A rT ° + r0' 2 + lu(cf))e x r (14) 
r 2 



1/ = - K*e x rTl + rcj)' 2 - \u{(j>)e x r 



a" = - (^-^ + -\<?'- uo 2 e v - x a - 2a{<f>)<f>W + 2A 2 (<P)e x W'(a 2 )a 



Here the components T$ and T^ are given correspondingly by 

T ° = ^ 2 A 2 (0)e"V 2 + l -A\4>)e^a' 2 + A\<P)W(a 2 ) (15) 



T\ = -\?A\$)e-»o 2 - l -A 2 ^)e~ x a' 2 + A\^)W(a 2 ) (16) 



The system (ITj) has to be solved at the following boundary conditions. We de- 



mand asymptotic flatness which means that z/(oo) = 0. On the other hand 
nonsingularity at the origin requires A(0) = 0. Concerning 0, nonsingularity at 
the origin implies <p'(0) = while at infinity <fi has to match the cosmological 
value 4>{oo) = (poo. Nonsingularity of a at the origin implies cr'(0) = O.We require 
finite mass an therefore we put o~(oo) = O.In addition we have to give the central 
value cr(0). 

It's well known that the tensor mass is just the ADM mass in Einstein frame 
(we note that in Einstein frame all mass definitions coincide). Therefore we can 
write directly the explicit expression for the tensor mass using the first equation 
of (|14D, namely 

Mt = 2g:L drr2 { K * T ° +e ~ x<p + 2 u((P) j = 2g;jo drr2v{r) (17) 

as T>(r) is define by the expression itself.Respectively, the particle number is given 
by 

N = An / drr 2 e^ (ujA 2 ((f))a 2 e-^) (18) 

When we are interested only in zeronodes solutions the boson star states are 
parameterized by the a field central value cr(0) or equivalently by p = a 2 (0). 
Let's consider two infinitesimally nearby configurations parameterized by <t(0) 
and a(0) + (5cr(0).The corresponding variations of the tensor mass and particle 
number are 



SM T = —— \ drr 2 5V (19) 

2G* jo 



SN = 4tt drr 2 5 (e*uA 2 {(j))cT 2 e-*) = (20) 

Air / drr 2 e*5 [ujA 2 ((j))a 2 e~^J + Air / drr 2 (uM 2 (0)o- 2 e _ 2 J 5e^ 
It's not difficult to obtain 



t^ / ^7^/r.2j 



$ e 2 = — e 2 x / drr 2 5V (21) 

2r Jo K ' 



Now substituting (pl|) in (^Cj) and after some algebra we have 

5N 



H = f°° d rr 2 e^8 LA 2 U)a 2 e^) + (22) 

An Jo v ' 

1 P OO /* oo 

+ -/ drr 2 5Vj d^(uoA 2 {4>)a 2 e-'i)ei x 
Taking into account the explicit form ([H]) of X> one obtains 

8 (uA 2 {<j>)a 2 e^) =^-U (^) + ^V* (^ 2 (0V 2 )) (23) 

-^ (iA 2 (0)e-V 2 + A 4 (0)WV) + le~ A 0' 2 + -^£7(0)) 
a; \2 K* 2/c* / 

In more detailed form the expression (^) is written as follows 

8 LA 2 ((f ) )e 2 e-t) = —5 (-) - — (lM + -e^'dA (24) 

— — 1 1 

-— (L a ba + A 2 ((P)e- x a'5a') - — ( -A 2 (0)a' 2 + —(f)' 2 ) 5e~ x 

CO V ' U) \2 K* J 

where L^ and L a are given by 

L+ = a{<f>) {-Lo 2 e- v A 2 {<\>)a 2 + A 2 ^y- X a' 2 + 4A 4 (0)W \a 2 )) + ^-U\<f>) (25) 

L CT = -a; 2 e-^ 2 (0)a + 2A\cf))W l \a 2 )a (26) 

Putting ( p4"|) in the first integral of (E^), performing integration by parts and 
taking into account the third and the fourth equation of the system ( |14|) one 
arrives at 

r°° \ / „\ r°° P2^ u +^) /7)\ 

/ drr 2 e^5 uA 2 (<P)a 2 e-^)= drr 2 - 5 — ) - (27) 

JO v J Jo U) Vft*/ 

_ / rfrr 2^ ( -A 2 (0)a' 2 + -0' 2 ) 5e" A 

JO Co" \2 K* J 

Substituting now ( P?D in (|2"2"| ) as one uses that 

£e~ A = — / drr 2 5£> 



r Jo 
we have 



- = l *r-A(r)*(-| 



where A(r) is the following expression 

A(r) = k* r dtt (±e^ + V ( l -A 2 {4>V 2 + -<f>' 2 ) + JuM 2 (0)<tVM x ) + (29) 



Using the first two equations of flli]) , it's not difficult one to show that A is 
actually a constant which turns out to be A = -.Therefore, we obtain finally 

= _ / drr 2 5— (30) 

4:71 UJ JO K* 



Comparing this expression with (19) we conclude that the following important 
relation holds 

5M T = u5N (31) 

The relation (|3l| ) may be rewritten in the form 

SM T 5N 

17 = u ¥ (32) 



Taking into account that < u) < m, it follows from (|32D that if p cr u is a critical 
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point for My (i.e. -M 1 - = 0), then p CT n is also critical point for N. 



Let My has a maximum at p = p crit ( s 2 T < 0).Then we obtain 
fS 2 N\ 1 (5 2 M T \ 



y„U ^l^„ <0 (33) 

which shows that A" has a maximum at p — p cr u, too. 

Therefore the tensor mass Mt and the particle number A" peak, as functions 
of the central density, at the same location . 

That the tensor mass and particle number peak at the same location results 
in a cusp in the bifurcation diagram Mt versus N. In infinitesimal small neigh- 
borhood of a cusp we may take expansions in powers of (p — p cr it) 

M T = M Tcrit + i i^0\ (p ~ Per*) 2 + \ {^0\ (P - Pcruf + 0(4) (34) 



1 {5 2 N\ , , 2 1 f6 3 N" 



A = N mt + ~ I ^ ) (P- Pcrit? + g I JT J (P - Pc"t) 3 + ( 4 ) 

\ < / crit \ " J cr it 



Now using (|32|) , the coefficients ( S M 2 T j and f g A f T j may be expressed by 
AT" = f ff ) and N'" = fff ) , and the result is 

\&P Jcrit \ 5 P 6 )criS 



'5 2 M^ 

crit 



2 - = ^critA^" 



Sp : 

7 




I -j-p- I = ^critA^ + 2u crit N . 
On the other hand, from the expansion for N we have 

p - Pcrit = ±((-^j(N crnt -N)j +\JL 
Substituting this expression in the expansion for Mt we obtain 

M T = M Tcrit + u crit (N - N crit ) t -^o'cnt fc) 2 (N mt - N)% + 0{2) (35) 

Let's denote by M^P and My correspondingly the tensor mass on the upper and 
lower branch of the curve Mt{N). Then making use of (J35|) one obtains 

M^ - M^ = --u/^t (^j 2 (iV crit - AT)I (36) 

It's easy to see that for the binding energy we have the same relation 

/critf^) 2 (Ncrit- N)i (37) 



9 /— 9\ I 

rpup t-iIow _ "I I * > 



These relations are scalar-tensor boson star versions of the similar relations in 
the fermion stars theory in pure general relativity [fL4| . It should be noted that 

3 - 

such dependence ~ (N crit — N)? is typical for catastrophe theory flTf . 

Conclusion 

Scalar-tensor theories of gravity violate the strong equivalence principle. This 
results in the appearance of three different possible masses as a measure of the 
total energy of the boson star. The stability analysis of the boson stars requires 
that the mass and particle number peak, as functions of the central density, at 
the same location. In this article we have proved analytically that it's the tensor 
mass which possess the desirable property. Therefore, it's the tensor mass which 
should be taken as the physical mass, for example in the construction of the 
binding energy of the star. While the numerical calculations have been done for 
boson stars in pure Brans-Dicke theory of gravity, our proof holds for the most 
general scalar-tensor theory.Especially, our proof involves a potential term for 
the gravitational scalar $. This is important, because many scalar-tensor models 
of gravity involve such term. On the other hand the potential term may play 
significant role in the early universe and to influence the boson star formation 
and stabilityalthough at present there aren't numerical investigations of boson 
stars in scalar-tensor theories with a potential term. 
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Finally, we believe that the result in this letter has a more general nature. It's 
the tensor mass which should be consider as the best candidate for the physical 
mass as a measure of the total energy of space-time in scalar-tensor theories of 
gravity. 
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